Two-dimensional Coulomb Systems in a disk 
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We consider two-dimensional Coulomb systems confined in a disk with 
ideal dielectric boundaries. In particular we study the two-component 
plasma in detail. When the coulombic coupling constant F = 2 the 
model is exactly solvable. We compute the grand-potential, densities 
and correlations. We show that the grand-potential has a universal 
logarithmic finite-size correction as predicted in previous works. This 
logarithmic finite-size correction is also checked in another solvable 
model: the one-component plasma. 
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1. INTRODUCTION 

Solvable models of Coulomb systems have attracted attention for quite 
some time. Very recently the two-dimensional two-component plasma, a 
model of two species of point-particles with opposite charges ±g at inverse 
temperature (3 = l/fc^T, has been solved in its whole range of stability F : = 
< 2 by using a mapping of this system into a sine-Gordon field theory. '0' 
With this mapping the grand-partition function and other bulk properties 
of the systern can be computed exactly. Mso some surface properties near a 
metallic wall's^ and an ideal dielectric wall® have been investigated. However 
this mapping into a sine-Gordon model does not give (yet) any information 
on correlation functions. 
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When the couphng constant F = 2 the corresponding sine-Gordon model 
is at its free fermion point and additional information on the system can 
be obtained. This fact is well known and much work has been done on the 
two-component plasma at F = 2 since the pioneer work of Gaudin on this 
model. In particular the two-component plasma at F = 2 has been studied 
in different geometries: near a plane hard wall and a dielectric interface^', 
a metallic wall®, in a disk with hard walls'^' and in a disk with metallic 
walls'^. . . . 

However it was not until very recently that the special case of ideal di- 
electric boundaries (that is Neumann boundary conditions imposed on the 
electric potential) has been studied by Jancovici and Samaj's-' for a system 
near an infinite plane wall or confined in a strip. The technical difficulty with 
this kind of boundary conditions is that the two-component plasma must be 
mapped into a four-component free Fermi field instead of a two-component 
free Fermi field as in all other cases of boundary conditions. 

A very interesting result of reference is that when the system is con- 
fined in a strip of width W made of ideal dielectric walls, the grand-potential 
per unit length exhibits a universal finite-size correction equal to tt/2AW 
which is the same finite-size correction for a system confined in a strip made 
of ideal conductor walls.® These finite-size corrections have been explained® 
by noting that if one disregards microscopic detail the grand-partition func- 
tion of a Coulomb system is the inverse of the partition function of the 
Gaussian model. 

Due to this analogy with the Gaussian field theory, a Coulomb system in 
a confined geometry is expected to exhibit logarithmic finite size-corrections, 
for instance in a disk of radius R the grand-potential should have a correction 
(1/6) In i?. For the analogy with the Gaussian field theory to be complete 
one should impose conformally invariant boundary conditions to the electric 
potential, for instance Dirichlet boundary conditions (ideal metallic walls) or 
Neumann boundary conditions (ideal dielectric walls). The case of Dirichlet 
boundary conditions was treated in reference 0. 

In this paper we study the two-component plasma at F = 2 in a disk with 
Neumann boundary conditions. One of the main motivations for this work is 
to show that the systems exhibits in fact the expected logarithmic finite-size 
correction. 

The outline of this paper is the following. In section 2, we present the 
model and adapting the method of reference |p we map the two-component 
plasma into a four- component free Fermi field. In section 3, we compute 
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the grand-potential and its large- i? expansion. We also compute the large- 
R expansion of the free energy of the one-component plasma which was 
solved some time ago.'S* In section 4, we compute densities and correlation 
functions. 



2. THE MODEL 



The system is composed of two species of point-particles with opposite 
charges ibg. The particles are confined in a disk D of radius R. It will be 
very useful to work with the complex coordinate z = re^'^ of a point r. The 
material outside the disk is supposed to be an ideal dielectric. This imposes 
Neumann boundary conditions on the electric potential: the interaction po- 
tential v{r, r') between two charges located at r and r' is the solution of 
Poisson equation 

Arv{r, r') = -2n6{v - r ) (2.1) 

with Neumann boundary conditions. However, it is well-known*0' that any 
solution of Poisson equation ( |2.1D in a closed domain D cannot satisfy ho- 
mogeneous Neumann boundary conditions dnv{r, r') = for r G dD, since 
Gauss theorem implies that ^gj^dnv{r,r') dl = — 27r. A natural choice is 
to impose inhomogeneous Neumann boundary conditions to the potential 
dnv{r,r') = — 27r/L, with L the perimeter of the domain D. In the case of 
the disk of radius R this gives 

dnv{r,r') = -1/R for r e dD (2.2) 

This impossibility for the electric potential between pairs to obey homoge- 
neous Neumann boundary conditions is not a problem for a system globally 
neutral in which the total electric potential will satisfy homogeneous Neu- 
mann boundary conditions if the pair potential satisfies ( p. 21) . It should be 
noted that for an infinite wall boundary it is possible for the electric potential 
between pairs to obey homogeneous Neumann boundary conditions.® 

The solution of Poisson equation ( |2.1| ) with boundary conditions ( |2.2| ) in 
a disk can be obtained obtained by usual methods of electrostatics (images, 
etc . . . ) . The solution is 

^(r, rO = - In '"~"^"f (2.3) 
where a is some irrelevant length scale. It can be easily checked that solu- 



tion (2^) can also be obtained as the limit of the boundary value problem 
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where outside the disk there is a dielectric with dielectric constant e — > (up 
to some constant terms).'"' 

It will be useful to introduce the image r* of a point r by z* = The 
electric potential between pairs can then be written as 



\z — z'\\z'* — z\\z'^ 



v{v, r') = - In J ^- ^ (2.4) 

and can be interpreted as the potential in r due to a charge in r' and a image 
charge with equal sign and magnitude located at r'*. 

If the temperature is not high enough the two-dimensional two-component 
plasma is not well defined, this is due to the collapse between pairs of op- 
posite sign. If the coupling constant T := [3q^ < 2 the system is stable. 
Here we will study the case F = 2. In order to avoid the collapse we will 
work initially with two interwoven lattices U and V . Positive particles with 
complex coordinates {ui} live in the sites of lattice U and negative parti- 
cles with complex coordinates {vi} live in the sites of V . We shall work in 
the grand-canonical ensemble and will only consider neutral configurations. 
The Boltzmann factor of a system with N positive particles and N negatives 
particles at F = 2 is 

N 

-I3Hm _ „2N 



X 



n(i?'-ki)(i?'-i^.i) (2.5) 



2 



The first product corresponds to the self-energies of the particles and the 
other terms to the interactions between pairs. Introducing the images this 
can be rewritten as 



J2 N N 



n(-^) IK-- (2-6) 



' 1=1 



Ui<jiUi - U,)iu* - U*)iUi - U*){U* - Uj){Vi - Vj){v* - V*){Vi - V*){V* - Vj) 
UijiUi - Vj){u* - V*){Ui - v*){u* - Vj) 



By using Cauchy's double alternant formula 



with 



(2.8) 



we find that the Boltzmann factor can be written as a 2N x 2N determinant 



-f3H„ 



UiVi 



det 



v 



1 



Ui - Vj Ui - Vj 



(2.9) 



Introducing the factors Ra/ui into the last rows of the determinant and 
the factors Ra/vi into the last columns, we finally arrive at the expression 



-ly det 



aR 



\ 



Ui - Vi 



aR 



UiVj — R^ 



(2.10) 



\ R^ — UiVj Vj — Ui j 



The grand-canonical partition function with position dependent fugacities 
C('u) and ^(f ) for positive and negative particles reads 



N 



7V=1 iil<---<«jv6(7 \j=l 
ui<---<ujySV' 



-(iHN 



(2.11) 



We shall now closely follow reference |^ to show that the grand-partition 
function can be written as a ratio of two Pfaffians. Using the same notations 
as in reference 0, let us introduce a couple of Grassmann anticommuting 
variables ('?/'^, 'j/'^) for each site m G f/ and similar Grassmann variables for 
each site in V . The Grassmann integral for any anti-symmetric matrix A 



de exp ( ^QM 



(2.12) 



with *0 - ' "''2 
Pfaffian of the matrix A 



. . . , . . .) and d9 = n. di^li^l Uu i'H^l is the 



= Pf(A) 

Let us denote the average of a quantity O by 



1 



dO O exp 
5 



(2.13) 



(2.14) 



It is well known tliatB 

{e,e,) = {A-% (2.15) 

For our purposes let us define the matrix A as having inverse elements 








uv 



aR \ 



u — V 
aR 



uv — R^ 
1 



\ i?2 - 

/ a 



uv V — u 



{A 



vu 



J 

aR \ 



V — u 
aR 



vu — R^ 

a 



(2.16a) 
(2.16b) 

(2.16c) 



(2.16d) 



vu u — V 



J 



The indexes [a, (3) G {1,2}^ label the rows and columns respectively. The 
matrix A is clearly anti-symmetric as required. We now introduce the anti- 
symmetric matrix M 



uu' 



vv' 















C(^) 





vu 

The Grassmann integral 

Z = J d9 exp 

is equal to 



Z = Pf (A + M) 
The ratio Z/Zq can be expanded in powers of the fugacities as 



N 



N 



^ = 1 + E E I nc(-oc(-.) ) ( n 



(2.17a) 

(2.17b) 

(2.17c) 
(2.17d) 

(2.18) 
(2.19) 

(2.20) 



Ar=l ui<. .<ujyG(7 \j=l 



li=l 



Using Wick theorem for fermions we find that 



N 



( 



n 



1)^ det 



11=1 



aR \ 



Ui — Vj UiVj — R^ 
aR a 



(2.21) 



\ R'^ — UiVj Vj — Ui J 



Comparing equations ( |2.20|) and ( p.21|) with equations (|2.11| ) and ( p.lO| ) we 
immediately conclude that the grand-canonical partition of the Coulomb 
system S is the ratio of two Pfaffians 



„ _ Z _ Pf (A + M) 
Pf (A) 



(2.22) 



Using the fact that the Pfaffian is the square root of the determinant we can 
write the grand-potential as 



(3n = -InS = -^lndet(l + K) = -^Trln(l + K) 
where the matrix K is K = MA^^ and has matrix elements 



(2.23) 



/ aR 
\ 



R^ - ^' 



— zz z — z 
a aR 



- zz' 



(2.24) 



We have introduced the notation u 
C(^) = C-(r) and (s, s') G {-, 



z — z 

(r,+), V = (r,-), Ciu) = C+(r) 



3. THE GRAND-POTENTIAL 



3.1. Formal expression of the grand-potential 



To compute explicitly the grand-potential from equation ( p.23| ) we must 



find the eigenvalues of K. From now on we will consider the continuum 
limit where the spacing of the lattices U and \^^oes to zero. In this limit 
it is natural to work with the re-scaled fugacity® m = 2TxaC/S where S is 
the area of a lattice cell. Also in this limit some bulk quantities will be 
divergent, because of the collapse of particle of opposite sign, and must be 
cutoff. Correlations in contrast have a well defined continuum limit. 
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Let {ipi"\r)}s=±-a=i,2 be the eigenvectors of m~-^K and 1/A the corre- 
sponding eigenvalues. The eigenvalue problem for K is the set of integral 
equations 



A 

2^ 



dr' 



D 



R 



1 



(2), 



— zz' 



dr' 



D 



z' — z 



z' — z 
R 



R^ - zz 



(2)/ 



(3.1a) 
(3.1b) 



These integral equations ( p.l|) can be transformed into differential equations 
plus some boundary conditions using the well-known identities 

d 1 d 1 



a - a- — = 7r5(r-r'). 

oz z — z' oz z — z' 

Applying dz to equation (|3.1a|) and dz to equation ( p.lbj ) yields 



2^ 



(2), 



dz^^Pir) 
d-z4'\r) 



(3.2) 

(3.3a) 
(3.3b) 



These differential equations ( p.3| ) can be combined into the Laplacian eigen- 
value problem 

A^/,(") = \^^p(<^) (3.4) 

which must be complemented with the following boundary conditions. If 
r = R is on the boundary, z = Re^^ , it can be easily seen from integral 
equations ( |3.1|) that 

^«(R) + e^Vf (R) = (3.5) 
An elementary solution of equation (p.4|) in the present disk geometry is 



with, from equation ( p.3b| ), 

V'S(r) 



(3.6) 



(3.7) 



where In is a modified Bessel functions of order £. Boundary conditions (|3.5| 
yield the following homogeneous linear system for the coefficients Ag 



- A^h+i{\R) + A+h{\R) = 
A.IeiXR) - A+I,+,{XR) = 



(3.8a) 
(3.8b) 
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For this system to have non trivial solutions its determinant must vanish. 
This gives the equation for the eigenvalue A 



From equation (|2.23|) the grand-potential then reads 



+00 



£=-oo 



+00 



£=0 



(3.9) 



(3.10) 



where the product runs for all A solution of equation ( p.9|) . The last equality 
in equation ( |3.1CI| ) comes from noticing that a change 1 leave 

equation ( p.9|) invariant. To evaluate the product in equation ( p.lO|) , let us 
introduce the analytic function for I positive 



f,{z) = {h{zRf-h^,{zRf) (^(^ 



2 

'zR 



(3.11) 



The zeros of this function are the eigenvalues A and it can be checked that 
f^{z)/ {zfe{z)) ^ as 2; ^ cxD, so this function can be factorized as a Weier- 
strass product 



(3.12) 



This function satisfies fi{0) = 1, /^(O) = 0, and fe{z) = fi{—z) so its zeros 
come in pairs of opposite sign and as a consequence the exponential factors 



in Weierstrass product (|3.12|) cancels out. We finally have 

= n(i - ^A) 



(3.13) 



where the product runs over all A solution of equation ( pl9| 
Then the grand-potential can finally be written as 



+00 



pn 



1=0 

+00 



mR 



2£ 



(£!)2 {l,{mRf - Ie+,{mR)') 



(3.14) 



The above expression is divergent and must be cutoff to a 



R/a where 



a is the ratio of the particles.*" 

It is interesting to notice that the grand-potential can be written as the 
sum of two terms 



ideal (m) + ideal {—IJl) 



where 



ideal (m) = -2^1n 



l=Q 



mR 



£\ {h{mR) + h+i{mR)) 



(3.15) 



(3.16) 



is the grand-potential of a two-component plasma at F = 2 confined in a disk 
with ideal conductor boundaries'i^ and 



ideal (-m) = -2 Vln 

r.on d ' " 



1=0 



mR 



l\ ihimR) - h+iimR)) 



(3.17) 



is formally the grand-potential of the two-component plasma with ideal con- 
ductor boundaries but with the sign of the fugacity changed (which of course 
does not correspond to any physical system). This interjif sting decomposition 
of the grand-potential also exist in the strip geometry.' 



■ (R) 



3.2. Finite-size corrections 

We now compute the large- i? expansion of the grand potential ( p. 14 ). 
First using decomposition (|3.15|) we can use the results of reference [§] for 
the expansion of the grand-potential with ideal conductor boundaries 

/3f}ideai(m) = -PpbTiR'^ + p-fc2nR+-\ii{mR) + 0(1) (3.18) 

cond Q 

where the bulk pressure is 

Pp,= !fL — + l) (3.19) 
ZTT \ ma J 

and the surface contribution with ideal conductor boundaries is 

f]^^ = ml -—In— ^ + -1 (3.20) 
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The second contribution to the grand-potential can be written as 



ideal (— m) 

2 cond 



where 



hard 



R/a 
R/a 

-2^ In 



1 - 



h+i{mR) 



h{mR) J 2 



h{mR) 



mi? 



(3.21) 



(3.22) 



is the grand-potential of a two-component plasma in a disk with hard wall 
boundaries. '-Q' The asymptotic expansion of this term was computed in refer- 
ence 171 and reads 



/3^^hard = -l3pbTiR^ + f3ih2TiR + - ln(mi?) + 0(1) 

o 

with the surface contribution for hard walls 

f^lh = m 

Finally, we only need to compute the asymptotic expansion of 

Ii+i{mRy 



1 1 
4 ~ 2^ 



(3.23) 



(3.24) 



S 



R/a 
i=0 



h{mR) 



(3.25) 



This can be done with Debye asymptotic expansions for the Bessel functions. 
First let us write 5* as 



R/a 



If{mR) niR 
The Debye asymptotic expansion for In/^ isAB 
In (mi?) 



(3.26) 



-iln(27r)-iln((mi?)2 + £2) 



+0 ((mRf + ey^i-^) 



mR 



(3.27) 
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Therefore 



h{mR) 



mR 



+ 



mR 



+ 



2 (mi?)2 + & ^{mRy + e mR 
imR 



e+^{mRy + 



(3.28) 



[(m/2)2 + £2]i/2 
+0((mi?)2 + f)-i/2) 

Using expansion (|3.28| ) and using Euler-MacLaurin formula for the sum over 



we finally find 



fix) rfx + -[/(0) + /(C,J] + 



S = mR +"^111 — + 0(1) 
2 mcr 



(3.29) 



(3.30) 



We notice that the divergent (when the cutoff a — » 0) surface contribution 
from (3 fl ideal (m) is canceled by the contribution from S giving for the surface 

tension the already known " finite expression 



m 



(3.31) 



Putting all terms together 



pn = -ppbTrR^ + p-fd27^R + - \n{mR) + 0(1) 



(3.32) 



The grand-potential has the expected (l/6)ln(mi?) finite-size correction. 



The one-component plasma 

As a complement to the above study of the finite-size corrections, in this 
section we consider another solvable model of Coulomb system, the two- 
dimensional one-component plasma at F = 2.^0 t!P This systems is composed 
of N particles with charge q living in a neutralizing uniform background. The 
one-component plasma in a disk with ideal dielectric boundaries was solved 
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by Smith/B The canonical partition function reads 



N 



N 



,3Ny4j^-N{N+l)/2 



(3.33) 



J] (7(/, N) - iV-(2^^+i-2')^(2iV + 1-1, N)) 



1=1 



with 'y{s,x) = t^~^e~^ dt the incomplete gamma function and Ath = 
h I \p2/KkBT is the thermal wavelength of the particles. 

We want to study the large-i? expansion of the free energy of the one- 
component plasma. The free energy can be written as 



f3F = (3F^,,^-J2^n 



n=0 



1 _ JY-(2iV-l-2n) 



-f{2N - n,n) 
-f{l + n,N) 



(3.34) 



where 



hard 



-3A^V4-iVln 



TxRa 



+ (A^(iV + l)/2)lnA^ (3.35) 



N 



-^ln7(n, A^) 



n=l 



is the free energy of a one-component plasma in a disk with hard walls 
boundaries.'"' The finite-size expansion of this terms is® 

/9Fhard = PfTrR" + /57 °- 27ri? + ^ In [{rrny^^R] + 0(1) (3.36) 
with the bulk free energy per unit "volume" (surface) 



_27r2a2 

and the "surface" (perimeter) contribution 



(3.37) 



n 




— / In 

27r Jo 



1 + erf (y) 



dy 



(3.38) 



n = N/{tiR^) is the density and erf(?/) = (2/y^) e dx is the error 
function. 
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The expansion of the remaining term in equation ( ^.34 ) can be obtained 
with the following uniform asymptotic expansions for the incomplete gamma 
function 



7(n + l,iV) 
7(2A^-n,A^) 



ni 
~2 



1 + erf (^^1^) +0[i/Vn' 



{2N-n-l] 



X 



1 + erf 



n- N 
V2N 



+ o (i/Vn) 



(3.39a) 



(3.39b) 



These expansions are valid when iV — n is of order y N and the corresponding 
terms in the sum (pj.34D are the ones that give a relevant contribution to (3F. 



Also for n such that — n is of order v N using Stirling formula we have 
the expansion for the factorials 



{2N -n-l)\ 



N\N 



N-n-l 



N\ N' 



n-N 



1 + fi/Viv' 



Finally replacing the sum in equation ( p.34| ) by an integral we find 

2erf(y) 



(3F = (3F^^,^ - ^2N I In ^ , „ , , 

/o l + erf(?/) 



dy + 0{l) 



(3.40a) 
(3.40b) 

(3.41) 



Putting this last result together with the expansion ( |3.36| ) for the hard wall 
case we find 

/?i^hard = /3/vri?2 + /57°-27ri? + ^ In [{T^nfl^R] + 0(1) (3.42) 



with 




= -v|^ / lnerf(y)rfy 



(3.43) 



For this model we find again the expected logarithmic finite-size correction. 



4. DENSITY AND CORRELATIONS 
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4.1. Green functions 

We return to the study of the two-component plasma. We are now inter- 
ested in the density and correlations functions. These can be obtained with 
the Green function 

as explained in reference The density ns{r) of particles of sign s is 

^s(r) = |$^Gr(r,r) (4.2) 

a 

and the truncated two-body density is 

nS(ri,r2) = E G:;.?(ri,r2)G--(r2,ri) (4.3) 



From its definition ( |4.1| ) the Green functions obey the integral equations 

G_,,(r,r2) = (4.4a) 



m 



R 


1 


R^-ziz 


Z — Zl 


1 


R 


Z — Z\ 


R^—ziz 


R 


1 


R'^-ziz 


Z — Zl 




R 



G-..(ri, r2) + — / rfr ( ^'^^^ '^-jf^ ) ^.(r, r2) 



z— Zl R^—ziz 
R 



1 

2^ 



R'^-ziZ2 Z2-ZI 
1 R 



(4.4b) 

Z2—Z1 R^ — ziZ2 / 

These integral equations can be transformed into the differential equations 
G-..(ri,r2)-^( ^° ^G,,(ri,r2) = (4.5a) 

)G_(ri,r2)-^a.(ri,r2) = -h{r, - r,)I (4.5b) 

where I is the 2x2 unit matrix. These equations can be combined into 

AriG'ss(ri, r2) - ^^^^^(ri, rg) = -m5(ri - ra)! (4.6) 



The boundary conditions can be obtained from the integral equations ( ^.4| ). If 
ri = R, ^1 = Re'^'^^ , is on the boundary then from the integral equations ( [4.4|) 
it can be seen that 

GL^(R,r2) + e^'^^GL^(R,r2) = (4.7a) 
GlUR,r2) + e'f^GfA^,V2) = (4.7b) 
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For the present disk geometry we look for a solution of equation ( [4.6| ) as a 
Fourier series in 0i. The solution for GH and G'^l can be written as 

Ztt ^ — ^ 

+AeIe{mri)Ie{mr2)] (4.8a) 

GL'(ri,r2) = ^5^e^^<^^5,J,(mrO (4.8b) 
Ztt ^ — ^ 

where Ki is a modified Bessel function of order i, r< = min(ri,r2), r> = 
max(ri,r2) and Ai and 5^ are constants (with respect to ri) of integra- 
tion that will be determined by the boundary conditions ( [4.7|) . From equa- 
tion ( |4.5a| ) we have 

G\Uri,r2) = ^d,,Gf,{rur2) (4.9a) 

GlUrur^) = -d,,Gll{r,,r,) (4.9b) 
m 



Therefore, using equations ( |4.8| ) we have, for ri > r2, 

Ztt ^ — ' 

+Aele-i{mn)] (4.10a) 

G'L'(ri,r2) = Ve^'*^5,_i/,_i(mri) (4.10b) 
Ztt ^ — ' 

Using the boundary conditions (O) we find the following linear system of 
equations for Ai and Bi 

AeIe{mr2)Ie + Be-iIe-i = -me-'''^'KeIe{mr2) (4.11a) 

AJ,{mr2)Ii.i + B,^J, = me-''^^K,_J,{mr2) (4.11b) 

where Ii = Ii{mR) and similar definitions for the other Bessel functions 
without argument. The solution of this linear system is 

A, = -me-^"f^ Keh + K,^iIi-i ^^^^^^ 



h{mr2) 
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And finally, 

Gllir^,r,) = iro(m|ri-r2|) (4.14a) 
Zn 

, K^I^ + Ki_ili_i ie(d,,-d,o) 
+ — Z^ 72 — —P2 h{mri)h{mr2)e"'^^^ 



G'}Jv^,r,) = ^^i:iiLKi(m|ri-r2|) (4.14d) 



27r |ri - r2 



—2 Ii^i{mri)Ii{mr2)e^ '^^ 



From equations ( [4 .41 ) it can be seen that the remaining Green functions can 
be easily deduced since they obey 



GL'.(ri,r2) = GLMri,r2) (4.15a) 
GLMri,r2) = GiMri,r2) (4.15b) 

Other useful symmetry relations between the Green functions are 



G'LMri,r2) = GiMr2,ri) GfAv.^v^) = Gf,{v2.v,) (4.15c) 
GL^(rl,r2) = -GL^(r2,rl) GLMri,r2) = -G^A^^.r^) (4.15d) 



4.2. Density 

The density is obtained from equation ( [4.2|) and it reads 

n,(r) =716 + — j2 /£(mr)' (4.16) 



where is the bulk density of the infinite system which is formally divergent 
when the cutoff vanishes. Writing formally the bulk density as 
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rearranging the terms in equation ( ^4.16| ) and using the Wronskian 



leKe-i + h^iKe = —- (4.18) 
mn 



the density at the boundary can be written as 

The above expression clearly vanishes since the term i is canceled by the 
term + So we recover the expected result 

UsiR) = (4.20) 

The strong repulsion between a charge and its image cause the density at 
the boundary to vanish. 

In reference |p it was shown that the density near an infinite ideal dielec- 
tric plane wall is 

2 

TTl 

ns{y) = nb-—Ko{2my) (4.21) 

/TT 

where y is the distance from the wall. To compare our result for the density 
in a disk and result (|4.21| ) near an infinite plane wall we plot in Figures 
and ^ both densities as a function of the distance from the wall for disks of 
different sizes R = 1/m and R = 3/m. In both cases the density for small 
distances y behaves as 

2 

ns{y)-nf, = —\n{2my) + 0{l) (4.22) 

ZTT 

It can also be seen in Figures |1| and ^ that the density decays faster for the 
semi-infinite system (plane wall) than in the disk case. This effect is stronger 
for the small disk mR = 1. 



4.3. Correlations 



From the Green functions ( |4.14| ) we obtain the two-body correlation func- 
tions using equation ( [4.3|) . Using the symmetry relations (|4.15| ) the correla- 
tion between a particle of sign s at ri and a particle of sign s' at r2 reads 



n 



(2)T/ 



iri,r2j 



-m 



G 



11 



ri,r2 



- G 



21 



ri,r2j 



(4.23) 
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Figure 1: Difference between the charge density and the bullc charge density 
ns{y) — as a function of the distance y = R — r from the wall. The dashed 
curve represents case of a two-component plasma near an infinite plane wall. 
The sohd curve represents the density in the disk case with mR — 1 



19 



-4 



Figure 2: Difference l:)etween the cfiarge density and the bulk charge density 
ns{y) — nfe as a function of the distance y = R — r from the wall. The dashed 
curve represents case of a two-component plasma near an infinite plane wall. 
The sohd curve represents the density in the disk case with mR — 3 
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This gives 



n 



ir(ri,r,) 



m 



Ko{m\Yi - ral) 



(4.24a) 



+ (2^1) 



E 



72-72 ^ 



and 



n>_is (ri,r2) 



4 

m 


^2 


- 2^1 


(27r)2 




- 



Ki{m\ri - rsl) 



(4.24b) 

2 



V2^j 



72 



/^(mri)/^(mr2) i^(0,_<^2) 



From equation ( [4.23| ) and the boundary conditions ( |4.7| ) it is clear that if one 
point is on the boundary 



n^^J, (ri, r2) =0 if ri G dD or r2 G dD 



(4.25) 



as expected due to the strong repulsion between a charge and its image. 
If r2 = the above expressions ( [4. 241 ) simplify to 



nSr(r,0) 



Ko[mr) H — Ii[mr) 



1 2 



(4.26a) 



+ 



(/2 - /2)2 



2\ 2 r 



2^7 



-KAmr) H -^5 — Iiimr) 



1 2 



(4.26b) 



/ m \2 /o(mr)2 

"12^^ (/2 - /2)2 
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Figure 3: Two-body density nL (r, 0) with one point fixed on the center of 
the disk for mR — 1. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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Figure 4: Two-body density nL (r, 0) with one point fixed on the center of 
the disk for mR = 3. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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Figure 5: Two-body density nljg (r, 0) with one point fixed on the center of 
the disk for mR = 1. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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Figure 6: Two-body density nijg (r, 0) with one point fixed on the center of 
the disk for mR = 3. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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It is interesting to compare these expressions with the bulk correlations for 
an infinite system'EP 

^SLkW = -(^^J[Mmr)f (4.27a) 

^SLik(r) = (^y[iri(mr)]^ (4.27b) 

Figures ^ and § show the two-body density n^^J'^{r, 0) for particles of same 
sign compared to the bulk values for different values of R and Figures H 
and y show the two-body density n_!,g (r, 0) for particles of different sign. 
For a small disk with mR = 1 there is a notable difference. In the disk case 
the correlations decay faster than in the bulk. This can be easily understood 
since there is a strong repulsion between a particle and the boundary. But 
this difference can be hardly noted if the disk is larger. For mR = 3 it can 
be seen in Figure § that the difference between the bulk and the disk case is 
very small (notice the change of scale in the vertical axis between Figures ^ 
and I). 

In Figures |^ and ^ we plot the structure function (charge-charge correla- 
tion) 

S(ri, r2) = 2{nfJ^{r,, r,) - ^^^(ri, r^)) (4.28) 

with one point fixed at center of the disk r2 = 0. For the small disk mR = 1 
there is a clear difference between the bulk case and the disk case. Due to 
the repulsion between a particle and its image, the screening cloud is more 
concentrated in the center of the disk than in the bulk case. But for the large 
disk mR = 3 the difference is hardly noticeable. Notice again the change of 
scale between Figures |^ and ||. It is interesting to note that there is not much 
difference between the bulk and the disk case if the radius of the disk is a 
few orders the screening length m~^ and larger. We only notice differences 
when R ~ m~^ and smaller. 

5. CONCLUSION 

We have studied the two-component plasma with coupling constant F = 2 
confined in a disk of radius R with ideal dielectric boundaries: the electric 
potential obeys Neumann boundary conditions. The model is solvable by 
using a mapping of the Coulomb system into a four-component free Fermi 
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Figure 7: Structure function S{r, 0) with one point fixed on the center of 
the disk for mR — 1. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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Figure 8: Structure function S{r, 0) with one point fixed on the center of 
the disk for mR — 3. The dashed curve represents the bulk correlation and 
the solid curve the correlation for the disk case. 
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field. We have computed the grand-potential, the density and correlation 
functions. 

The grand-potential can be formally written as an average of the grand- 
potential for ideal conductor boundaries and the same grand-potential for 
ideal conductor boundaries but with the sign of the fugacity changed. For 
ideal conductor boundaries the surface tension is infinite when the cutoff 
vanishes. Here, the average makes the surface tension finite. This fact also 
appears for a two-component plasma in a strip, 

The Neumann boundary conditions for the electric potential being confor- 
mally invariant it is expected that the grand-potential of the system exhibits 
a universal finite-size correction (1/6) Ini?. This was explicitly checked on 
this solvable model. We also checked this universal finite-size correction on 
the model of the one-component plasma ^t F = 2 in the same confined ge- 
ometry which was solved some time ago.^tJ 

The density vanishes for a point on the boundary of the disk. This is 
expected since there is a strong repulsion between the particles and the 
boundary due to the image forces. This is also true for the correlations, 
they vanish if one point is on the boundary. We compared the correlations 
functions for an infinite system without boundaries and the present system 
in a disk. Due to the repulsion between the particles and the boundary, the 
screening cloud around a particle in the center of the disk is smaller than the 
one for an infinite system. But the difference between the screening clouds 
is only noticeable for disks with radius of the order of the screening length 
and smaller. If the disk has a radius a few orders larger than the screening 
length, the difference can be hardly noted. 
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